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High-energy emission from the pulsar striped wind: a 
synchrotron model for gamma-ray pulsars 
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ABSTRACT 

Gamma-ray pulsars constitute a class of high and very high-energy emitters for which 
the known population is steadily increasing thanks to the Fermi/Large Area Telescope. 
Nowadays, more than a hundred such pulsars have been detected, offering a reasonable 
sample onto which to apply statistical techniques in order to outline relevant trends in 
the averaged properties of this (maybe not so) special class of pulsars. In this paper, 
their gamma-ray luminosity and spectral features are explained in the framework of 
synchrotron radiation from particles located in the stripe of the pulsar wind. Apart 
from radiative losses, particles are also subject to a constant re- acceleration and re- 
heating for instance by a magnetic reconnection induced electric field. The high-energy 
luminosity scales as L 7 « 2 x 10 26 W(L sd /10 28 W) 1 / 2 (P/l s)" 1 / 2 where L sd is the 
pulsar spindown luminosity and P its period. From this relation, we derive important 
parameters of pulsar magnetosphere and wind theories. Indeed, we find bulk Lorentz 
factor of the wind scaling as T v » 10 Tree (^sd/10 28 W) 1 / 2 , pair multiplicity k related 
to the magnetization parameter a by K G Tree 

« 10 8 , and efficiency r\ of spin-down 
luminosity conversion into particle kinetic energy according to the relation r\ a ks 1 . A 
good guess for the associated reconnection rate is then r rec ~ 0.5 (L s d/10 28 W)~ 5 / 12 . 
Finally, pulses in gamma-rays are visible only if L s d/P > 10 W/s. This model dif- 
fers from other high-energy emission mechanisms because it makes allowance not only 
for rotational kinetic energy release but also for an additional reservoir of energy an- 
chored to the magnetic field of the stripe and released for instance by some magnetic 
reconnection processes. 

Key words: Acceleration of particles - magnetic fields - radiation mechanisms: non- 
thermal - MHD - gamma-rays: theory - pulsars: general 



1 INTRODUCTION and comparison between radio and gamma-ray radiation if 

both are available. 

Since their discovery in the radio band more than 40 

years ago, pulsars have now been firmly detected in all It was long suspected that polar cap or outer/slot 

wavelengths of the electromagnetic spectrum, from radio gaps could explain this emission. Nevertheless, recent Fermi 

through optical, X-ray and eventually high-energy gamma- obser vations clearly d isfavored the polar cap explana- 

rays. Recently the Cr ab was even detected in very-high en - tion jAbdo et alj|2010h . Some of these gamma-ray pulsars 

ergies by VE RITAS (IVERITAS Collaboration et al.l l201lh do not significantly show a spectral cut-off around a few 

and MAGIC (|Aleksic et alj|201lh . Although the radio emis- GeV but rather a significant change in the spectral index 

sion mechanism remains poorly understood, the MeV- of the power law. Moreover MAGIC/ VERITAS detection 

GeV pulsed emission puts severe constraints on the high- of the Crab above 100 GeV seems to rule out outer gap 

energy counterpart radiation models. The new catalog of models because the accelerating electric field combined to 

gamma-ray pulsar s obtained by the Fermi-LAT instrument radiation reaction limited flow renders it difficult to observe 

|Abdo et afll2010h increased the number of gamma-ray pul- photons above a few GeV. The determination of the precise 

sars from seven to about fifty. Since then, new pulsars are location of the emission regions is still problematic. The in 

disco vered regularly, m ore than hundred are listed nowa- phase pulsation observed i n both gamma-ray and in radio for 

days (|Nolan et al]|2012f ). This allows for the first time a rea- some millisecond pulsars (|Guillemot et al.ll2012h furnishes 

sonable statistical analysis of the high-energy emission prop- one more apparent contradiction between radio and high- 

erties of these objects like spectral shapes, cut-off energies, energy emission mechanisms. Do they eventually originate 
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in the same place in the magnetosphere or in the wind? But 
why are then some gamma-r ay pulsars radio quiet as those 
seen by iPletsch et all l|2012l ). We are still far from a con- 
sensus on the emission model, be it in the radio band or in 
gamma-rays. 

Radio pulses and gamma-ray photons are expected to 
be produced in different emission sites, probably close to the 
neutro n star surface for the former, des cribed by a polar cap 
model (Radhakrish nan fc Cooke 1969), and in the vicinity 
of the light-cylinde r for the latter, explained by outer gaps 
l|Cheng et al.lll986l ). 

Recently, gamma-ray light-curves have been com- 
puted for a realistic magnetospheric model based on 
3D MHD simulations o f the near pulsar magnetosphere 
l|Bai fc Spitkovskvl l2010). In this model, gamma rays are 
expected close to the light-cylinder. 

An alternative site for the production of p ulsed radi- 
ation has been investigated a few years ago by iKirk et all 
(2002). This model is based on the s triped pulsar wind, 
originally introduced bv lCoronitil (|l990l ). Emission from the 
striped wind originates outside the light cylinder and rela- 
tivistic beaming effect s are responsible for the phase coher- 
ence of this radiation (Petri 20091 l201ll ) 

In this paper, we show that the pulsed high-energy emis- 
sion up to a few GeV as observed by Fermi/LAT can be 
explained by synchrotron radiation emanating from the rel- 
ativistically hot current sheet present in the pu lsar striped 
wind. Our approach follows the study done by iLvubarskiil 
(1996) who included possible pair creation which are dis- 
carded in the present work. Sec. [2] describes the dynamics 
and geometry of the relativistic plasma flow as well as the 
radiation model. Then Sec. [3] derives some important con- 
straints on the bulk Lorentz factor of the wind, the pair 
multiplicity factor, the magnetization parameter as well as 
a criterion for gamma-ray pulsar detection. Conclusions are 
drawn in Sec. [4] 



2 THE SYNCHROTRON STRIPED WIND 
MODEL 

Synchrotron emission from the striped wind applies suc- 
cessfully to t he pulsed optical po larization properties of the 
Crab pulsar jPetri fc Kirkll2005l ). Our aim is to extend to 
the highest possible energies the synchro-photons emanat- 
ing from the current sheet. In the present work, we are not 
concerned with the phase-resolved emission nor in the polar- 
ization properties, but we only focus on the average gamma- 
ray features of the Fermi/LAT detected pulsars depending 
on two fundamental observables: the period P of the pul- 
sar and its first derivative P. We emphasize that this two- 
parameter family of model remains very restrictive and will 
not be able to explain in detail each individual pulsar. In- 
deed, the mechanisms occurring in the current sheet are far 
from being understood, as they involve micro-physics that 
has not yet been addressed self-consistently and linked to the 
overall structure and dynamics of the magnetosphere and 
wind. Nevertheless, as we will show, our model is able to re- 
produce qualitatively and quantitatively the full set of Fermi 
data with reasonable accuracy. The derived parameters such 
as wind Lorentz factor, pair multiplicity and magnetization 
are consistent with values obtained from other independent 



considerations, giving us confidence in our description of the 
emission mechanism. 



2.1 Plasma configuration 

At the heart of any pulsar emission model, be it radio, opti- 
cal or X-rays/gamma-rays, there is a (possibly relativistic) 
plasma flowing in a strongly magnetized field, preferentially 
within the magnetosphere for the polar cap and slot/outer 
gaps. Contrary to these models, we assume that the pulsed 
emission is produced within the current sheet of the striped 
wind, thus well outside the magnetosphere, a place where the 
magnetic configuration switches from poloidal dominant to 
toroidal dominant. The plasma in this current sheet is the 
crucial ingredient in our model. It is embedded in a mostly 
non emitting cold plasma. More precisely, the striped wind 
flow is made of two distinct parts, namely 

• a cold and strongly magnetized plasma of particle den- 
sity number n c and Lorentz factor r v , as measured in the 
lab frame. Note that F v corresponds to the bulk Lorentz 
factor of the whole striped wind structure. This part of the 
wind does not radiate significantly. 

• a hot and tenuous but weakly magnetized plasma of 
particle density number rij, occupying a fraction A of the 
wavelength of the wind A v = 2 n fi v tl , also measured in 
the lab frame. The fraction of hot plasma is negligible with 
respect to the cold part, i.e. A -C 1. Tl = c/Q is the light- 
cylinder radius, /3 V = v/c the wind speed normalized to the 
speed of light c and Q the stellar rotation rate. 

Later on, it will be useful to deal with the proper densi- 
ties and proper lengths given by n c — T v n' c , — T v n' h 
and At = T v A v . Proper quantities in the rest frame of the 
wind are always primed except for thermodynamical quanti- 
ties such as pressure and temperature which are meaningful 
only in the proper frame. Furthermore, we assume that pairs 
are created within the magnetosphere and cool down quickly 
before reaching the wind zone. Therefore, they exclusively 
replenish the cold part of the striped wind. Acceleration and 
heating of particles occurs once in the stripe and not be- 
fore. This is clearly an additional source of energy, a kind 
of magnetic luminosity, independent of the spin-down lumi- 
nosity of the pulsar. Therefore, rotational kinetic energy is 
spent to produce relativistic leptons with a curvature radia- 
tion reaction limited Lorentz factor of about 7 ~ 10 7 . They 
will cool down in the magnetized part of the wind because 
of magnetic field strength around the light-cylinder higher 
than 10~ 2 T acco rding to Fermi/LAT first pulsar catalog 
l|Abdo et al.l |2010| ). Indeed, the synchrotron cooling time 
T syn ~ 7.7 S7 _1 (y^p) is much less than the period of 
a pulsar. In the most unfavored case, radiation losses are 
still such that r syn < 1 ms. Thus, our model can apply irre- 
spective of the period of the pulsar, millisecond or normal. 



2.2 Wind flow 

Although the hot plasma is much less populated than the 
cool part, it radiates significantly more because of its rel- 
ativistic temperature. Indeed, to ensure a quasi-stationary 
equilibrium state, magnetic pressure outside the stripe has 
to be compensated by gaseous pressure inside the current 
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(i) 



sheet. Thus in the wind frame, pressure balance implies 

1 / / 2 B' 2 
- 7 h n h m c c = - — 
3 2 )i 

m c is the electron rest mass and uo vacuum permeability. 
We recall that the magnetic field is essentially toroidal in 
the wind zone, a Lorentz boost thus gives B' = B/T v where 
B — BlVl/v, Bl being the magnetic field strength at the 
light-cylinder. The Lorentz factor 7^ can be deduced from 
considerations about the dynamics in the current sheet as 
explained in the next paragraph (see Eq.©), leading to the 
knowledge of the density of the hot plasma n' h as given by 
the pressure balance Eq. (0. 

The bulk kinetic energy of the wind is extracted from 
the spin-down luminosity of the pulsar. Noting r\ the effi- 
ciency coefficient of rotational to particle energy conversion, 
the plasma kinetic energy density is 



F v n c m c c = 77 



47T r 2 c 



(2) 



Note that we neglect the contribution of the hot component 
because A <C 1 and L s d = An 2 I P P~ 3 is the spin-down 
luminosity and I = 10 38 kgm 2 the neutron star moment of 
inertia. Moreover, the cold plasma density is related to the 
pair creation rate N± (particles generated per second) by 
4 7r r 2 c n c = N± where the contribution from both poles are 
taken into account and 

3 



N± w 2.77 x 10 30 s" 1 k 



P_ 



( B n 



\ 10 s T / V 10 km 



Rn 



(3) 

k is the usual pair multiplicity factor, e the absolute value 
of the charge of an electron, eo the vacuum permittivity, 
B na the stellar surface magnetic field intensity and i? ns the 
neutron star radius. Inserting Eq. © into Eq. ©, we found 
a relation between F v , k, and r\ such that 



r vK Ri 8.7 x it) 8 ?? 



( L s d 

V 10 28 w 



1/2 



(4) 



A similar relation holds for the magnetization defined by 
a = B 2 /(fj,o F v n c m c c 2 ). The magnetization parameter is 
approximately the inverse of the efficiency coefficient, i.e. 

(1 + sin 2 x) °"n = 1 (5) 

where \ is the pulsar obl iquity and assumi ng the force-free 
spi n-down lum inosity of ISpitkovskvl l|2006l ) and confirmed 
bv IPetril Uml). This reduces the number of independent 
parameters in the model. As expected, the fraction of spin- 
down luminosity converted into particles remains weak r\ -C 
1 as long as the magnetization parameter remains large a 3> 
1. 



Particles in the sheet lose energy due to synchrotron 
cooling. However, a stationary state in the wind frame can 
be reached if particles re-energize by acceleration in a recon- 
necting magnetic field configuration. The key parameter is 
the collisionless reconnection rate, r rec . It will be constrained 
from spectral features of the Fermi/LAT data. This means 
that the induced electric field re-accelerating particles is of 
the order E lec « T Tec cB. Equaling the synchrotron cooling 
timescale with the acceleration timescale in the reconnection 
layer, we find that the maximum Lorentz factor of the lep- 
tons (we use primed quantity to express the fact that they 
are evaluated in the wind frame) is 



3 fioe c r 



2 a T B' L r L 



(6) 



Basically, it depends on the magnetic field strength at the 
light cylinder and on the distance to the centre of the pulsar. 
For these Lorentz factors, the typical energy of a synchrotron 
photons as seen in the wind frame is 



/ 3 ,2 B 
£s =2 7h 



2 _ 9 fi em c 
4 or B„ 



(7) 



where B q ~ 4.4 x 10 T is the quantum magnetic field. Thus 
the synchro-photon energy does only depend linearly on the 
reconnection rate. Interestingly the photon lies in the MeV 
band because e' B = 236 MeVr roc . Back to the lab frame, the 
photon is boosted by a Lorentz factor at most 2 F v thus 



e B = 2 T v e' B = 472 MeV T v r r< 



(8) 



According to the Fermi/LAT catalog (|Abdo et al hold ), the 
cut-off energy lies in the range 0.2-8 GeV. In our model, 
we interpret this cut-off more properly as a break in the 
spectral energy distribution. To make the idea more precise, 
think about an inverse Compton spectrum evolving from 
the Thomson regime to the Klein-Nishina regime. There is 
a break in the power law entailing precious information on 
the underlying emitting plasma. Here, the break is a conse- 
quence of the rapid cooling of electrons more energetic than 
7^ m c c 2 . Using a fiducial value of r rcc w 0.1, we see that the 
Lorentz factor of the wind should be around T v G [5, 10 2 ]. 
We emphasize that the estimates given here should be un- 
derstood as orders of magnitude roughly within a factor 10. 
More confidently we expect T v € [1, 10 3 ] as will be shown 
from the observations. Moreover, the reliability of the value 
of the cut-off energy is probably doubtful at least for the 
most energetic pulsars like the Crab if not for almost the 
whole sample of gamma-ray pulsars. 



2.3 Dynamics in the current sheet 

It is usually claimed that the total gamma-ray luminosity 
should not exceed the spin-down luminosity. Actually, this 
assertion only holds if the rotation of the neutron star is the 
unique source of energy. Another non negligible reservoir 
could be the magnetic field itself. We will indeed assume 
that part of the radiation is due to magnetic reconnection 
within the current sheet. This will alleviate the restriction 
L 7 < L s< j although the model not necessarily violates this 
condition. It will be checked a posteriori. 



2.4 Synchrotron luminosity 

In the frame of the current sheet, relativistic particles in 
the hot component radiate synchrotron power at a rate 
P^ndLB') = (4/3)a TC1 £U' B where U' B = B' 2 /(2 Mo ) 
is the magnetic energy density in the wind frame. The 
total emissivity of the hot plasma, assuming isotropy is 
jsyn.(v') = {n' h I '(4 71-)) P s ' yn (7h, B') 5 (v' - u' syn ) where u' syn = 
(3/2) 7^ 2 e B' /m e . By Lorentz transformation and using rela- 
tivistic invariants, the total synchrotron luminosity, boosted 
to the observer frame (expected to be interpreted as the 
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gamma-ray luminosity) is 

T o t in- 20 ,„ /Tree A r 5/4 Lad / n s 

L syn = 8.7x10 WW— — -pLJt si— (9) 
v i v 

where T v > 1 is the bulk Lorentz factor of the wind, A < 1 
the thickness of the current sheet in fraction of the wave- 
length A v , and a > 1 sets the minimum distance where 
emission starts (in the region r > oa)- From eq. Q, we 
get 

L syn = 3.9 x 10" 7 W ^V /2 A (10) 

Knowing the synchrotron luminosity Eq. (|10p and the re- 
lations and iJSJl, we are able to deduce important con- 
straints on the fundamental parameters. Theoretically, for 
fixed {a, A, «, 77, r rcc }, our model predicts L syn oc y' L s d/P. 



3 RESULTS 

In this paragraph, we discuss the implication of the striped 
synchrotron model and the latest Fermi/LAT detections and 
give explicit values for the most relevant parameters of pul- 
sar magnetosphere and wind theories. 

3.1 Luminosity and spectral features 

It is possible to relate F v to the pair multiplicity k and to 
the magneti zation a. Accordin g to the second source cata- 
log (2FGL) jNolan et alj|2012h . a good fit of the whole sam- 
ple of gamma-ray pulsars following the formal dependence 
as suggested by Eq. (|10|) thus with respect to \J L s d/P is 

^ - 2 X 1026 W (l^) V2 (^)" /2 (") 

The constant of proportionality comes from Fermi/LAT 
data and will be compared to the factor in front of yj L s d/P 
in Eq. ()10|) . Note that the integral energy flux above 
100 MeV from 2FGL and distances fro m ATNF are sub- 
ject to several sources of uncertainties llManchester et all 
2005 J]. This reflects into the derived quantities such as 
{r v , Tree, A, a} for which the same remarks hold. We also 
used a beaming factor of 1. Thanks to the above fit, Eq. (|11[1 . 
and to the relation ()10|) . we get (n/rff/ 2 A y/r rcc /a « 
5 x 10 18 . Adopting typical values of A w 0.1, a w 10, 
the relation between pair multiplicity and efficiency becomes 
KTrJc ~ 1.2 x 10 8 77. From this, we immediately deduce the 
bulk Lorentz factor of the wind 

r ^ io ^ 5 (ib^r (i2) 

It is proportional to the square root of the spin-down lumi- 
nosity. By comparison between the cut-off energies observed 
by Fermi/LAT and the prediction Eq. p). we can constrain 
the reconnection rate for each pulsar such that 

/ 4.72 GeV y 5/6 / L sd y 5/12 
Trcc ~ \E cut (GeV) J \W™Wj 

1 see also the website from ATNF 

http: / / www.atnf.csiro.au / research / pulsar / psreat / expert.html 
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Figure 1. Pulsar gamma-ray luminosity L 7 versus spin-down 
luminosity L s( j. Fermi observations are depicted by red points 
whereas the predicted luminosity is presented by blue tri- 
angles. The magenta line corresponds to the law L 7 = 
10 26 W(L sd /10 26 W) 1 / 2 , the green line to L 7 = L sd and the 
cyan line to the best fit Eq. {TTJ . 

The reconnection rate is about unity for low L s d pulsars 
such as the millisecond ones. This might be possible if the 
background pair densit y is about 1% of the current sheet 
particle density number ()Bessho fc Bhattacharieef2012l ) . We 
suspect magnetic dissipation to reach its maximal possible 
value for those pulsars with L s( j < 10 2S W whereas r rcc de- 
creases monotonically down to 0.01 for the brightest pulsars 
with L s d > 10 28 W. This splits the gamma-ray pulsars pop- 
ulation into two groups, a first with reconnection pushed to 
its highest limit and a second with less drastic reconnec- 
tion rate. Moreover, the photon spectral indexes founded by 
Fermi/LAT fall in the range F £ [0.6,2]. This is consistent 
with a synchrotron spectrum emanating from a power law 
distribution of leptons with p G [1, 3] and dnh/d'y( r y) oc 7~ p . 
This could be explained by a relativistic kinetic reconnec- 
tion process, either very localized close to the X-point where 
it is found that p ~ 1, or on a global scale, within an ex- 
tended region such as a simulation box for PIC simulations 
ijjaroschek et al.ll2004l ) where p ~ 3. A detailed analysis of 
the kinetic rec onnection processes i s out o f the scope of this 
letter but see IZenitani fc Hoshinol |2007l ). Detailed knowl- 
edge of the particle distribution functions within the cur- 
rent sheet await s extensive PIC simulations with radiation 
effects included (|jaroschek fc Hoshino| [2009l. The results of 
our model concerning the expectation about the gamma-ray 
luminosity are summarized in the plot Fig. [T] We show the 
gamma-ray luminosity versus the spin-down luminosity for 
a large sample of gamma-ray pulsars. The predicted lumi- 
nosity, blue triangles, is overlapped to the Fermi/LAT data, 
red points. Although not explicitly imposed by our model 
and actually not required, the gamma-ray luminosity never- 
theless satisfies L 7 < L s d. 

3.2 Condition for pulsed gamma-ray emission 

Gamma-ray emitters are identified as pulsars because of 
their very stable clock. In our model, pulsation is a con- 
sequence of the combination of relativistic beaming and 
striped wind structure. Therefore, we expect to observe 
MeV-GeV pulsed emission only if the stripe survives the de- 
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Extinction line 
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Figure 2. The extinction line for gamma-ray pulsation in the 
PP diagram. All radio pulsars are shown in grey points whereas 
Fermi/LAT detected gamma-ray pulsars are in red triangles. 
They all satisfy our criterion Eq. 114> , lying above the blue ex- 
tinction line. 

structive influence of magnetic dissipation due to reconnec- 
tion. Assuming as a first approximation that the stripe re- 
sembles a relativistic Harris current sheet, a stationary equi- 
librium can exist only if the drift speed of the leptons within 
the sheet is less than the speed of light . This is actually a 
criteri on for reconnection as explained in lKirk fc Skiaeraasenl 
(2003). In the wind frame, the drift speed is given by j3' d = 
r' B /(A X' v ) where r' B is the proper Larmor radius. Replacing 
the bulk Lorentz factor of the wind by the approximate re- 
lation Eq. (|12J1 , we derive a simple criterion for pulsed emis- 
sion from the condition /3' d < 1. According to its spin-down 
luminosity and to its period, a pulsar will be detected as 
gamma-ray emitter if L Bd /P > (10 23 W/s) r rcc a 3 /A 2 . For 
typical pulsar parameters, we found 

^■ET ~ 1027 W /s ( 14 ) 

This criterion is shown in Fig. [2] for the Fermi/LAT pulsar 
second catalog and marked as a blue line. All detected pul- 
sars fall into this criterion. The two populations of pulsars, 
millisecond and normal, are clearly distinguishable. The gap 
between the extinction line and the detected millisecond pul- 
sar is due to the P _1 law in Eq. (|14p . The average parame- 
ters (a, A, r rec ) may also be slightly different of those from 
normal pulsars. The absence of pulsation does not imply 
that the underlying neutron star is not a gamma-ray emit- 
ter, it just says that it will not be seen as a gamma-ray 
pulsar. 



and bulk Lorentz factor of the wind T v are explicit func- 
tions of L s d. Moreover the predicted extinction line is in 
agreement with gamma-ray pulsar observations. 

Our synchrotron model predicts a clear break in the 
spectra around a few GeV depending on the bulk Lorentz 
factor and reconnection rate. In order to explain recent ob- 
servations by VERITAS and MAGIC, this picture should be 
supplemented with an inverse Compton counterpart. The 
model would also benefit from phase-resolved high-energy 
polarization measurements which are highly discriminating 
for current models. 
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4 CONCLUSION 

In this paper, we showed that the pulsed MeV-GeV emission 
from gamma-ray pulsars can be explained by a synchrotron 
model in the striped pulsar wind including some magnetic 
reconnection. Relevant average parameters of our model are 
given by A ~ 0.1 and a ~ 10. The reconnection rate r roc 
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